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We have investigated the potential energy surface for neutral and charged alanine 
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and glycine chains consisting of up to 6 amino acids. For these molecules we have cal- 

culated potential energy surfaces as a function of the Ramachandran angles tp and i/j. 
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Our calculations are performed within ab initio theoretical framework based on the 

density functional theory and also within semi-empirical model approaches. We have 
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demonstrated that the excessive positive charge of the system influences strongly its 

geometrical and conformational properties. With increasing of the excessive charge 
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amino acid chains become unstable and decay into two or more fragments. We have 

ON . 

analysed how the secondary structure of polypeptide chains influences the formation 

of the potential energy landscapes. We have calculated the energy barriers for tran- 
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sitions between different molecular conformations and determined the ones being 
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energetically the most favourable. 
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I. INTRODUCTION 

This work was presented on "The eighth European Conference on Atomic and Molecular 
Physics" (ECAMPVIII) (Rennes, France, July 6-10, 2004) an d on the "Electronic Structure 
Simulations of Nano structures" workshop (ESSN2004) (Jyvdskyld Finland, June 18-21, 

2004). 
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Amino acids are building blocks for proteins. Recently, it became possible to study 
experimentally fragments of proteinSji.e. chains of amino acids, in a gas phase with the use 
of the MALDI mass spectroscopy [llyLyl. Ab initio theoretical investigations of amino acid 
chains began also only recently J4J, y, la |7|, |s|, |3, |l(| and are still in its infancy. 

We have investigated the potential energy surface for neutral and charged alanine and 
glycine chains consisting of up to 6 amino acids. For these molecules we have calculated 
potential energy surfaces as a function of the Ramachandran angles ip and tjj often used for 
the characterization of the polypeptide chains |llLll2J|. Our calculations are performed within 
ab initio theoretical framework based on the density functional theory and also within semi- 
empirical model approaches. We have demonstrated that the excessive positive charge of the 
system influences strongly its geometrical and conformational properties. With increasing 
of the excessive charge amino acid chains become unstable and decay into two or more 
fragments. We have analysed how the secondary structure of polypeptide chains influences 
the formation of the potential energy landscapes. We have calculated the energy barriers 
for transitions between different molecular conformations and determined the ones being 
energetically the most favourable. 

II. THEORETICAL METHODS 

Our exploration of the potential energy surface of alanine and glycine chains is based on 



the density-functional theory (DP 
semiempirical AMI method 



/_(DFT]ac 



accounting for all electrons in the system and on the 



Within the DFT one has to solve the Kohn-Sham equations, which read as (see e.g. 

ft \ 

y + U ions + V H + V xc J ipi = Eiipi, (1) 

where the first term represents the kinetic energy of the z-th electron, and Ui ons describes its 
attraction to the ions in the cluster, Vh is the Hartree part of the interelectronic interaction: 



Vn(rl= ljf^-dr>, (2) 



and p(r') is the electron density: 



N 



where 14c is the local exchange-correlation potential, ^ are the electronic orbitals and N is 
the number of electrons in the cluster. 

The exchange-correlation potential is defined as the functional derivative of the exchange- 
correlation energy functional: 

_ SE xc [p] 
Vxc ~ Sp(r) ' [V 

The approximate functionals employed by DFT methods partition the exchange- 
correlation energy into two parts, referred to as exchange and correlation parts. Both parts 
are the functionals of the electron density, which can be of two distinct types: either lo- 
cal functional depending on only the electron density p or gradient-corrected functionals 
depending on both p and its gradient, Vp. In literature, there is a variety of exchange 
correlation functionals. In our work we use the Becke's three parameter gradient-corrected 
exchange functional with the gradient-corrected correlation functional of Lee, Yang and Parr 
(B3LYP) lialiatLsl- We utilize the standard 6-311++G(d,p) and 6-31G(2d,p) basis set to 
expand the electronic orbitals ipi. 

III. RESULTS OF CALCULATION 

In figure^we show the dihedral angles tp and i)j that are used to characterize the potential 
surface of the polypeptide chain. 

In figure |21 we present the optimized geometries of the alanine and glycine polypeptide 
chains that have been used for the exploration of the potential energy surface. All geometries 
were optimized with the B3LYP density functional. We used the 6-31++G(d,p) and 6- 
31G(2d,p) basis sets to expand the electronic orbitals in the molecule. 

In figures I3TIT21 we present the potential energy surfaces for the polypeptide chains pre- 
sented in figure 121 



In figure EI] we show the optimized structures of the alanine tripeptide. Different geome- 
tries correspond to the minima on the potential energy surface (see contour plot in figure 

El). 

In figure ITU we show the optimized structures of the glycine tripeptide. Different geome- 
tries correspond to the minima on the potential energy surface (see contour plot in figure 

EJ). 

In figure El we show the optimized structures of the glycine hexapeptide in helix confor- 
mation. Different geometries correspond to the minima on the potential energy surface (see 
contour plot in figure EJ) . 

The geometries of singly charged alanine and glycine dipeptides are shown in figure fTHl 
The ionization of the system changes dramatically its potential energy surface and the 
secondary structure as it is seen from the contour plots presented in figures El and El 



In figure El we show the distribution of observed dihedral angles tp, ip of non-G 



ycine 



residues in protein structures selected from the Brookhaven Protein Data Bank 11 
The circles show conformations corresponding to the forbidden regions. 
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FIG. 1: Dihedral angles tp and ip that are used to characterize the potential surface of the polypep- 
tide chain 
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FIG. 2: Optimized geometries of the alanine and glycine polypeptide chains calculated with the 
B3LYP/6-31++G(d,p) (a and b) and B3LYP/6-31G(2d,p) (c-f). a) Alanine tripeptide; b) Glycine 
tripeptide; c) Alanine hexapeptide in sheet conformation; d) Glycine hexapeptide in sheet confor- 
mation; e) Alanine hexapeptide in helix conformation; f) Glycine hexapeptide in helix conforma- 
tion; The energies below each image are given in a.u. 
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FIG. 3: Potential energy surface for the glycine tripeptide calculated with the B3LYP/6-31G(2d,p) 
method. Energies are given in eV. Numbers mark energy minima on the potential energy surface. 
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FIG. 4: Potential energy surface for the alanine tripeptide calculated with the B3LYP/6-31G(2d,p) 
method. Energies are given in eV. Numbers mark energy minima on the potential energy surface. 



6xGlycine (Helix); B3LYP/6-31G(2p,d) 
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FIG. 5: Potential energy surface for the glycine hexapeptide in helix conformation calculated with 
the B3LYP/6-31G(2d,p) method. Energies are given in eV. Numbers mark energy minima on the 
potential energy surface. 
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6xGlycine (Helix); AMI 
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FIG. 6: Potential energy surface for the glycine hexapeptide in helix conformation calculated with 
the AMI method. Energies are given in eV. 
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6xGlycine (Sheet); B3LYP/6-31G(2p,d) 
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FIG. 7: Potential energy surface for the glycine hexapeptide in sheet conformation calculated with 
the B3LYP/6-31G(2d,p) method. Energies are given in eV. 
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6xGlycine (Sheet); AMI 
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FIG. 8: Potential energy surface for the glycine hexapeptide in sheet conformation calculated with 
the AMI method. Energies are given in eV. 
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6x Alanine (Helix); B3LYP/6-31G(2p,d) 
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FIG. 9: Potential energy surface for the alanine hexapeptide in helix conformation calculated with 
the B3LYP/6-31G(2d,p) method. Energies are given in eV. 
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6*AIanirie (Helix); AMI 
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FIG. 10: Potential energy surface for the alanine hexapeptide in helix conformation calculated 
with the AMI method. Energies are given in eV. 
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6xAhinine (Sheet); B3LYP/6-31G(2p,d) 
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FIG. 11: Potential energy surface for the alanine hexapeptide in sheet conformation calculated 
with the B3LYP/6-31G(2d,p) method. Energies are given in eV. 



16 



6xAlnnine (Sheet); AMI 
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FIG. 12: Potential energy surface for the alanine hexapeptide in sheet conformation calculated 
with the AMI method. Energies are given in eV. 



17 



1 

E =-818.4465 (a.u.) 




'^=-156.629 

'F= 163.643 




<P=-166.121 
»l'=-54.8826 




AE= 0.2296 (eV) 



AE=0.3424 (eV) 



c 

<t =^78.31 9 
1'= -51 .1003 



<P =-82.7052 
* = 75.6576 



A E =0.2928 (eV) 




ip=69.7371 
T=-31.3 



6 

AE=0 2992 (eV) 



,*>*£> 



<P =59 8445 

>l'=37.8723 



<P = 83.2Z29 
H' =-53.1411 



8 

AE=0.3834 {eV) 



v^*V 






(('=598552 
'l'=-133.143 



FIG. 13: Optimized structures of the alanine tripeptide. Different geometries correspond to the 
minima on the potential energy surface (see contour plot in figure^). 
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FIG. 14: Optimized structures of the glycine tripeptide. Different geometries correspond to the 
minima on the potential energy surface (see contour plot in figure [3]). 
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FIG. 15: Optimized structures of the glycine hexapeptide in helix conformation. Different geome- 
tries correspond to the minima on the potential energy surface (see contour plot in figure EJ) 
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FIG. 16: Optimized geometries of singly charged alanine (a) and glycine (b) tripeptides calculated 
with the B3LYP/6-31G(2d,p) method. Energies below each image are given in a.u. 
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3 x Alanine, Charge +1 ; AMI 
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FIG. 17: Potential energy surface for the singly charged alanine tripeptide calculated with the 
AMI method. Energies are given in eV. 
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3xGlycine, Charge +1; B3LYF/6-31G<2p,d) 
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FIG. 18: Potential energy surface for the singly charged glycine tripeptide calculated with the 
B3LYP/6-31G(2d,p) method. Energies are given in eV. 
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FIG. 19: Distribution of observed dihedral angles <£ , ^ of non- Glycine residues in protein structures 
selected from the Brookhaven Protein Data Bank 
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